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The  marginal  success  of  l  ho  several  h  j  gh-reso.)  u  l  i  on 
frequency-wavenumber  (-f-k)  techniques  to  data  is  cited  from 
the  literature.  Their  ability  to  resolve  signals  from  two 
closely  spaced  sources  is  not  markedly  superior  to  that  of 
ordinary  beamforming.  Moreover,  such  nonlinear  techniques 
yield  distorted  magnitudes  and  azimuths.  The  ordinary  -f- k  ' 
^spectrum^  is  shown  to  be  no  more  than  a  1-signal  estimator, 
and  the  high  resolution  techniques  to  be  but  variations  of 
that  1-signal  estimator.  In  this  paper  th<-  notion  of.  the 
wavenumber  spectrum'  is  set;  aside.  ny  analogy  to  the  1- 
signal  estimator  (the  ordinary  f-k  spectrum)  a  linear  M  - 
signal  estimator  is  developed.  The  high  resolving  power  of 
this  technique  and  the  fidelity  of  its  estimates  is  demon¬ 
strated  theoretically  and  by  computer  examples  both  real  and 


'V- 


synthetic . 


The  frequency-wavenumber  spectrum,  which  is  a  multi¬ 
dimensional  equivalent  of  the  ordinary  frequency  spectrum, 
is  used  in  the  sciences  for  theoretical  and  experimental 
analysis  of  traveling  waves.  It  was  introduced  formally 
into  seismology  by  Burg  (5)  in  an  appl  i.eation  to  data 
analysis.  The  ordinary  unsmoothed  three  dimensional 
frequency-wavenumber  spectrum  of  time  series  data  sampled 
at  discrete  points  in  space  is  given  by 


M  ft- 
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where 


is  the  j  iulo;:  o 
sample  points. 


f.  the  spacial 
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is  the  iinite  Fourier  trans¬ 
form  of  t he  nth  time  series. 

is  the  vector  wavenumber, 

1  s  the  vector  location  of  the  p  £ 
sonso!  ,  o i  sum p  1.  e  p o i  n  t  . 


Each  Fourier  transform  term  is  equivalent  to  n  sinusoid. 

For  example,  the  sinusoid  for  the  A,  t  k  transform  at  frequency 
<ju)  has  amplitude  AaC^)  and  phase 


(at  the  center  of  the  time  window) . 


Now  —  k.  * 


is  the  phase  delay,  between  the  origin 


and  f'n  ,  of  a  plane  wave  arriving  from  the  ;i?  i.muth  of 
the  vector  k  and  travel  ing  at  Lite  phase  velocity 


So,  multiplication  of  the  fransfuim  b\  the  kotnoL  C?  X.  p  ( I  k- 
has  the  effect,  of  advanc i ng  the  sinusoid  hy  'just  the  amount 
the  wave  itself  had  delayed  it.  Tlius  Lhe  summation  in  (1), 
above,  is  a  beam  sum,  and  the  I-  -  It  spectrum  is  just  the 
frequency  domain  equivalent  of  ordinary  beam  steering. 

When  the  traveling-wave  delays  are  exactly  compensated 
for  by  the  beam  shifting,  i.c.,  when  the  true  !\  of  the 
signal  is  selected,  the  sinusoids  add  up  constructively  with 
no  interference,  and  the  power.  P,  is  maximized.  Within 
certain  limits,  then,  maxima  or  peaks  .in  k  space  are 
treated  as  indications  of  the  presence  of  traveling  plane 
waves,  and  the  location  and  size  of  the  maxima  are  taken 
as  estimates  of  the  speed,  bearing,  .frequency  and  power 
of  those  signals.  If  more  than  one  signal  is  present  or 
if  there  is  noise  in  the  data,  though  exact,  determinations 
are  no  longer  possible,  the  Pk  spectrum  is  still  useful 
for  detecting  and  estimating  signals,  again  within  limits. 

One  of  those  limitations  is  .imposed  by  the  finite  width 
of  the  maxima  associated  with  sign:.  Is.  (  !)  )  The  case  is 
analogous  to  that  of  the  ordinary  i requency  spectrum  in 
which  components  are  represented  by  peaks  ot  finite  width. 


spectrum  have  a  half- 


Plane  wave 
power  w id tli 

* 
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signal  peaks  in  the 
of  the  order 


h-k 


AX 


where  AX  is  the  width,  or  aperature,  of  the  array  of 
spatial  sample  points.  (4)  If  two  signals  in  the  same 
time  vindow  and  frequency  band  are  also  close  enough  in 
phase  velocity  and  azimuth  so  their  wavenumbers,  say  Vv. ( 
and  kx  ,  are  such: 

k,-|g|  <  &i<, 

then  their  maxima  in  the  -f-k  spectrum  are  merged  and 
form  a  single  peak.  (23)  Thus,  because  the  sensor  arrays 
are  spatially  finite  their  resolving  power  is  finite.  Attempts 
to  increase  that  power  of  resolution  through  data  processing 
technique  have  required  mathematical  schemes  to  reduce  the 
width  of  the  lobe  of  the  signal  peak  (1-3  6-1. ‘p  17,  19). 

However,  the  straight-forward  geometric  appeal  of  this 
approach  has  proved  misleading  thus  far.  In  such  hybrid 
spectra  signal  lobe-widths  indeed  have  been  narrowed 
substantially.  Nevertheless,  when  signal  pairs  approach 
each  other  in  the  k-plane,  resolution  still  fails  as  the 
separation  nears  A  k,  to  wit  ,  the  lobe  halt: -width  for 
the  ordinary  -  k  spectrum.  (2,  11,  13,  19,  20). 
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Observations  of  other  i  nvcs  l  nju  Loi  s  on  t:)io  shortcomings  of 
various  nigh-resolution  f requency-vuvonuml >or  techniques  are 
cited  below. 

Lintz  (1968)  finds  chat  1  ho  r  ; -wv  r-  so  !  \u  i  on  -f  -  k  spectral 
technique  of  Haney  (1967)  does  not  .-...1  >mi  t  1  cant.  1  y  improve  the 
capability  of  a  seismic  array  to  detect  multiple  time-overlap- 
pint  events  from  different  azimuths. 

Salat  and  Sax  (1969)  expo rimon Lully  1 >nd  the  high-resol¬ 
ution  -f- k  spectrum  of  Haney  (19o7),  and  that  of  Capon  (1968), 
(1969)  no  better  at  resolving  two  simultaneously  arriving 
waves  than  the  ordinary  {-k  spectrum.  McGowan  and  Lintz 
(1968)  call  attention  to  an  unrecoverable  distortion  of  the 
true  amplitude  spectrum  in  Haney's  technique,  and  the  marked 
disadvantage  of  spurious  peaks  under  certain  conditions  which 
they  regard  as  the  inevitable  result  of  using  a  high-gain 
procedure . 

Seligson  (1970)  describes  conditions  under  which  Capon's 
high-resolution  technique  displays  less  "angular  resolution" 
than  ordinary  beamforming.  McDonough  (1972)  concludes  that 
variations  in  amplitude  from  sensor  to  sensor  may  be  expected 
to  produce  anomolous  behavior  in  Capon's  processor.  Of  course, 
just  such  variation  in  amplitude  from  sensor  to  sensor  will 
result  precisely  because  of  the  presence  of  two  or  more  signals 


McDonough  offers  arguments  to  show  that  ordinary  bsamforming 
is  less  susceptible  to  instability  resulting  from  small  signal 
modeling  errors  than  all  other  array  processors. 

Haney,  too,  notes  that  in  the  processor  he  describes 
(1967)  variation  in  amplitude  from  sensor- to-sensor  could 
distort  the  spectrum  beyond  recognition.  He  remedies  this 
difficulty  by  forcing  the  same  amplitude  upon  each  input 
channel,  thus  destroying  the  very  amplitude  information  that 
would  be  indicative  of  the  presence  of  two  or  more  signals. 

Woods  (1973)  concludes  that  given  favorable  conditions, 
the  resolving  power  of  the  maximum-likelihood  -f“k  spectrum 
can  be  effectively  infinite,  but,  disappointingly,  offers  com¬ 
puter  examples  on  synthetic  data  in  which  the  input  signal 
pairs  are  well  spaced  to  begin  with  (they  are  separated  by  a 
distance  of  0.9  of  the  main-lobe  half-width).  Cox  (1973)  also 
offers  theory  suggesting  that  given  arbitrarily  high  signal- 
to-noise  ratios  arbitrarily  fine  resolution  should  be  possible, 
but  he  does  not  offer  a  method. 

It  may  be  argued  that  the  limited  resolving  power  of  the 
several  high-resolution  techniques  results  from  the  wavenumber 
spectrum  being  in  reality  a  1-signal  estimator.  Indeed,  the 
ordinary  -f-k  "  spectrum"  is  a  least  squares  estimator  for 
fitting  given  data  to  a  single  plane  wave.-,  as  shown  further 


on.  In  routine  automated  processing  of  the  LAS A  LP  data 
Mack  and  Smart  (1973)  found  the  ordinary  spectrum  useful 
for  estimating  only  one  signal  at  a  time.  Estimates  of  a 
possible  second  signal  were  made  by  recomputing  the  wave- 
number  spectrum  after  the  first  (and  larger)  estimate  had 
been  subtracted  from  the  data.  They  call  this  process  strip¬ 
ping?  it  is  useful,  of  course,  only  for  estimating  signals 
separated  by  about  the  reciprocal  of  the  array  diameter  or 
more.  At  that,  such  estimates  of  a  pair  of  signals  are  not 
optimum,  but  first  order  approximations. 

Properly,  the  f-k  spectrum  is  defined  only  for  signals 
of  infinite  spacial  extent  traversing  infinitely  large  arrays. 

The  effect  of  a  signal  of  wavenumber  kB  is  then  confined 

- - 

to  the  point  k*  in  the  spectrum.  Approximations  to  this 
definition  are  useful  if  the  dimensions  of  signals  and  arrays 
are  sufficiently  large.  Failing  that,  the  "spectrum"  reduces 
to  a  1-signal  estimator  as  noted.  While  the  high-resolution 
techniques  do  attempt  to  extend  the  effective  array  diameter, 
they  all  test  the  wavenumber  space  with  a  ] -signal  probe,  as 
in  the  ordinary  -f- k  spectrum. 

It  is  proposed  here  to  set  aside  the  notion  of  a  spec¬ 
trum.  Rather  we  will  extend  the  1-signal  estimator  to  an 
M  -signal  estimator  thus  to  permit  the  simultaneous  removal 


of  the  effects  of  one  signal  from  the  estimate  of  another 
and  so  achieve  true  high-resolution.  t  the  same  time,  use  of 
beamf ormincj  (in  the  k  -plane)  to  estimate  each  of  the  AA 
signals  will  preserve  the  stability  and  estimate  fidelity  of 
the  ordinary  {-k  spectrum. 

In  the  following  discussion  a  1-signal  least  squares 
estimator  is  developed  and  is  identified  v/ith  the  ordinary 
|-k  spectrum.  Analogy  to  the  1-signal  estimator  is  used  to 
develop  an  M -signal  estimator. 


Conventional  Frequency-Wavenumber  Analysis 

In  the  conventional  frequency-wavenumber  spectrum 
(ordinary  or  high-resolution)  a  single  plane  wave  is 
hypothesized  at  each  frequency.  That  model  is  then  tested 
over  the  wavenumber  space  of  interest.  One  attempts  to 
minimize  the  error 

A/ 

<£■  -  22  iiv-Ae’"' 

r\a  / 

by  varying  A  and  k where 

are  the  complex  Fourier  series  terms 
(for  the  given  frequency) 

is  the  sensor,  or  channel,  index 

is  the  total  number  of  sensors 

are  the  location  vectors  of  the 
sensors 

is  the  complex  Fourier  series  term 
for  the  hypothesized  plane  wave 
(at  the  given  frequency) 

is  the  wavenumber  of  the  hypothetical 
plane  wave  (at  that  given  frequency) 
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r\  =  }  N 


is  the 


model,  i.e.,  the  hypothesized  plane  wave. 

Note  that  also  one  can  write  £ 


*  -  z:  u. 


_n  i  on 

-//<  V 


-A 


s  ince 
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For  a  given  is  minimized  by  setting  A  to 

a/  — -n* 


/A  ~  ‘a. 


a/  kJ  ^  Cr 

A«/ 
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which  is  shown  by  the  following: 

Let  J  CV\ 


=  UA  e-'7< 


<3  -/■  t  c  s 


Then 


<£■  = 


ZI  +i(Ca- c) 


ns/ 


=  z:  c 


-  <$)  +  C  c*  “  c) 


Take  partial  derivatives: 

A / 


£€■  ; 


eP<~ 


-af  (c -c) 


The  second  term  is  the  ordinary  frequency-wavenumber 
spectrum  yy  _ * 

Pc -fy  k)  ~  7s/  2T  Ujc/0 <^~l ^  ^ 


becomes  exactly  zero  when 


that  is,  when  the  data  describe  a  single  plane  wave  exactly. 
The  smaller  is,  in  a  given  situation,  the 

more  likely  is  the  hypothetical  plane  wave 

A 

because  the  smaller  drC  k.)  is,  the  larger  the  F-statistic 
is  for  the  hypothesis.  The  F-statistic  is  given  by 

F  =  (M-0-  P('^4c"k) 

This  single  plane  wave  model  is  often  applied  in 
attempts  to  analyze  a  2-signal  case  (or  a  possible  2-signal 
case) .  In  such  an  analysis  each  signal  is  treated  as  if  it 
existed  by  itself,  the  presence  of  the  other  being  ignored 
with  consequent  distortion  of  estimates  by  mutual  interfer¬ 
ence.  This  interference  can  be  serious,  and  if  the  two 
signals  are  not  separated  in  k.  -space  by  at  least  the  half¬ 
width  of  the  main  lobe  of  the  array  response,  they  are  likely 
to  appear  as  but  one  signal,  their  main  lobes  having  coalesced 
Attempts  to  improve  the  performance  of  the  single  wave  hypo¬ 
thesis  (in  application  to  the  two  signal  case)  have  been 
made  in  which  the  main  lobe  of  the  array  response  has  been 
slenderized  mathematically  by  alternative  methods  of  csti- 


mation  of  the  wavenumber  spectrum.  The  object  has  been  to 
reduce  the  main-lobe  half-width  and  so  resolve  signal  pairs 
which  otherwise  have  coalesced  main-lobes  indistinguishable 
from  a  signal  case.  These  results  have  been  marginal.  In 
the  various  high-resolution  techniques  the  influence  of  the 
one  signal  on  the  analysis  of  the  other  has  been  ignored, 

Analysis  of  possible  2-signal  cases  calls  for  a 
2-signal  model,  in  particular  when  the  2-signals  are  known 
(or  suspected)  to  be  so  close  together  as  to  have  their  main 
lobes  merged. 

As  the  1-signal  model  serves  for  both  the  0-  and  the 
1-signal  case,  so  one  might  expect  a  2-signal  model  to  be 
effective  in  all  three  cases:  0,  1,  or  2-signals. 


Multiple  Signal  Frequency-Wavenumber  Analysis 


By  analogy  to  the  1-signal  model,  one  would 


expect  to  solve  a  2-signal  model  by  minimizing  the  error 


A/ 


/w 


=  23 


varying  A,  kj  B,  and  'k  ,  where 


B 

$ 


is  the  complex  Fourier  series 
term  for  the  second  hypothesized 
plane  wave  (at  the  same  given 
frequency) 

is  the  wavenumber  of  the  hypothe 
tical  plane  wave  (at  that  same 
given  frequency) 


There 


are 


now 


two  signals  to  solve  for: 
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Again, 

let 

A  ~  c9  + 

/  c 

,  A*-  c3-/C 

Taking  first  partial  derivatives  while  noting  that 
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Setting 
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is  in  the  1-s.ignal  case, 
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Therefore, 
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Analogous ly 


B  =  TrZ:  (u„-Ae  k-r:) •  e- *  ^ 


In  this  form  A  and  &  are  optimized,  that  is,  they  produce 
thf  minimum  value  of  <*=  for  any  arbitrary  pair  of  k  and 
/l<  •  Adopting  the  notation: 


P s  -jg 22  L)„e"' 
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one  may  write  simply: 
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-(B*i|u./‘5  +bA^u:j^) 

+  (  A* A  +  B*E>) 

+( A*B  -kt^X)?'  +  ®i(l: 


^  =  Nf  £  U*u,  -  (  A*P4  AP*)  -  (  B',;a  4  BQ*) 
+  (  A*A  4  B*B)  +  (  A*  B  E  4  A  8*  E*) 

Rearranging  the  terms  in  , 

^  =  7j£  U*U,  -  (AA  +  AP*)  -(6*0  4  6  0.*) 

4  A*(  A  4  BE)  4  B” ( A E* 4  B) 

and  recalling  that 

P=  A  4  BE  and  Q  =  AE’-PQ 

^=Tr£,U*lA  - (AP*  +  BQ*) 


Further  bstituting 

A  »  P-  6E  and 


B  U  -  A  E 1 


x- 


^  — 


“  -(p*p  +  ctg  -  pq‘e'-  p'aE) 

n-i  N  " - / 
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or,  wut^en  out, 


*\=-\ 


The  identity  of  these  last  2  equations  may  be  demonstrated 
by  noting  that  the  numerator  (above)  equals 


=  P*P  +  Q*Q  - PQ* E*  —  P^Q E 


Since  ^  is  a  sum  of  squares,  by  definition  it  must 


be  non-negative  everywhere.  Therefore  the  second  of  the 
2  terms  in  ^  ,  above,  must  always  be 


-I-Tv 
N  /- 


Thus  to  minimize 


,  one  must  maximize 


A*‘  1  j  =  i  J  -]  - 1 

I  -  1 

A  =  I 

This  is  the  2-signal  test,  analogous  to  the  ordinary  frequency- 
wavenumber  spec trum,  which  is  the  1-signal  test.  However,  it 
is  more  convenient  to  retain  the  form 

AP*  +  BQ* 

This  2-signal  “f*“*  k  "spectrum”  then  is  computed  from  3 
beams  fas  the  ordinary  f-k  spectrum  is  computed  from  1  beam) . 
The  beams  are 


i  1X1 

=  n-£  U. 

'  r  s» 


-i  k-fc. 


i  lN1 

"FT 
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the  mean  of  the  data  transforms 
that  'nave  been  beamed  to  U. 

(one  of  the  two  wavenumber 
variables ) , 

the  mean  of  the  data  transforms 
after  beaming  to  £  (the 
other  wavenumber  variable), 


E_  _L  /t(k-kK£  Which  is  the  (complex)  array 

N  ^  response 

a=i 

This  2-signal  test  is  solved  as  is  the  ordinary  -f-k 
spectrum,  numerically,  by  searching  the  wavenumber  space  of 


interest.  Now,  however,  there  are  4  dimensions  to  search, 


over  which  to  test  the  error  criterion. 

It  is  instructive  to  submit  a  known  pair  of  pure, 
noiseless  signals  to  the  2-3ignal  test  to  illustrate  the 


function  of  the  elements  of  the  expression: 

Let  U  =  F  F  +  G 


)  *\  *  1  > 


Beaming  them  exactly  to  and  k  (since  these  are  known 
in  this  special  case), 

P  =  irZZ  (Fe1  k-‘R;  4-  G  P'1" N  •  “-1  P* 
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=  FE*  +  G 


Then 
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(  F+QE  -  (F£*  +  Q)E)/(l-el‘E) 

(  F  +  QE  -  FE*E-QE.)/0 -£*e) 
F(i  -  e*e.)/(i- e:*e)  =  F 


B=(FE4+G  -(F  +  GE)Ej/(l-E*E 


B  =  <3(i -e,‘e.)//(i-e^e)  =  G 


Thus 


I 

-  -kZ 


i 


-Be*' 


I 

-*£  f 


Q  e‘^'^  „  Fe^^-CJe'^  ' 


<£*  =  O 


This  little  exercise  clarifies  a  bit  the  function  of  the  array 
response,  El,  in  the  signal  models  A  and  B 

The  development  of  the  2-signal  test,  of  course, 
suggests  the  derivation  of  a  3-signal  test,  by  analogy: 

First,  the  form  of  the  test  would  be,  analogously, 

<s  -  TrZ  I U*  -  A  -  B  Jk*  -  C  2 


Introducing  the  notation 


R  —  n  UA  e 


“  N  e 


itk-k)*rA 


E,-TtZ=irt'k’’?  E> 


N  4 


and  expanding  e  : 


^  -  Tf|  (u*  -  A*-^  -  B':'o  -  C'&Aj 


*(ua-a  -  b 


«5r  ** 


-( AMfirZ:  +  A  -kt  Uje? l<"') 

a*  i 

's<  *\<  t 

+  C-kZU*^-) 
+  (A*A  +  B*B  +  c*c) 


+  (A"B^Z  e'(4- k)T"  +  A  B*-^'  e 1  ('k'-^R) 


Now  noting  that  in  the  2-signal  test 


1  ? 


so  that 


P  «  A  4  BE  and 

Q=  AE*+  E> 


P  E 
IQ  1 
I  E 
E*  I 


E>  = 


i  p 

E*  Q 

I  E 
E*  I 


one  recognizes  that,  in  the  3-signal  test, 


=  A  4  B  E,  4  CE, 

=  ae;‘+  b  +  ce* 

=  A  E"  +  B  Ej*  +  C 


and,  defining 


I  E,  E, 

=  E,  I  Ea 

E»  E»  I 


P  E,  E, 


A  -  9 


E.  •  d 


R  Ef  I 


e  A 


etc . ,  or 


A  =  [  P(  I  -E‘Ea)  4  Q(E*EI-E.)  4  R(E,Es-E$/de* 

B  -  t  P(e;e,-ED  +  Qd-E^EJ  +  RCE^-E^doA 

C  =  [  P(EfE;-E>“)  +  Q(E,’E?-6)  +  R(l-ETE)]^Ln 

ders=  I  -  E*E,  -  E*E,  -  EJEi  4  E,  E;'lE,  +  E*E,E« 


Now  rearranging 


N  LlflV  (a**P+  B^Q+CR)  -  (a  P*4  BCt*  4  CR*) 

n»t 

+  A*(A4  BE,  4  CE^  +  ^(AEf+  B  +  CE3)  4c’!(AE?+  B  Ej  ■*  C  ) 

and  substituting  P  ,Q  ,  and  R 

.  NJ 

•=^  =  7?^  U?LL  -(apm4  bq^cr*) 


To  minimize  ,  then,  one  will  maximize 


A  P*  +  BQ*+  C  R* 


the  3-signal  test,  or  3-signal  analog  to  the  conventional, 
1-signal  frequency-wavenumber  spectrum.  The  function  is 
composed  of  6  beams:  P  ,Q  ,  and  R  ,  the  3  beams  of  the 
data,  U  and  E.  ,  Ea  ,  and  Ej  ,  the  3  beams  of  the  array 


response , 


Remembering  the  1-signal  test  (conventions 


i  f-k 


spectrum) , 


i  l<  •  |  a 


*  =TTZ:U,-Ae’ 


=  Tli  (u*-AV’k'^(u^-Ae‘H 

A  -  t 

we  may  rewrite  it  as 

=  Tl£u.rU*-A*P-AP*  -  A*  A 


(since  P  *  ) 

A»l  ' 


or 


M 

■k  *  =  TiZ  U*U«  ~  A*A  -  A  P*+A*A 

A*l 
M 


-  wzuru,  -AP- 


A  =  l 


Thus 


ap*  •  *(*|£uk.nH‘) 


is  the  expression  one  must  maximize  in  order  to  minimize 
the  error.  So  the  f-k  spectrum  (for  the  1-signal,  conven¬ 
tional,  case)  is 


AP 


and 

and 


AP*  +  BQ.*  is  the  2-signal  test, 

AP*  •+  BQ*  +  C  R*  is  the  3-signal  test. 


In  the  1-signal  test 

A- 


For  the  2-signal  test 

Ip  e,| 


A- 


a 


I  E, 
Ef  I 


For  the  3-siynal  test 


A  - 


B  = 


P  E,  E* 
Q  l  E, 
R  Ef  I 


I  P 

E,rt  Q 


IE, 
Ef  I 


I  E  E 
ETI  E 
E^E?I 


e.i  c. 


*3M«.7nU’S5* . r ia; V.^^f 


fe<  ft 


a  & 


l\\ 

Hi 

I'l 

I ! 


This  formalism  makes  evident  the  relationship  between  the 
successive  tests.  Thus  one  may  extropolate  and  directly 
write  the  expression  for  the  M-signal  test  in  simple,  terse 
form.  For  example,  the  4-signal  test  is 

AP*+  BQVCR*+DS* 

in  which  S,  the  sum  of  the  data  beamed  to  yet  a  4th  point, 
^  ,  is  introduced  into  the  sequence,  P,  ,  and  P  ; 
and  in  which 


A  = 


P  E,E*E. 
Q  I  E5E5 

r  e;i  e. 
s  e:e: 


ll 

.1  _  •* 

1  { 

I'l  :■ 

p  ■ 

;« 

and  t~~  ,  ^ 

|  ’ 

(f-k) 

IEEE 
E,  I  E,E 

e;e:  i  e 
e:  e:  e: 


1  etc . 


is  the  array  response  at  <&-lo .  E  5  ,  that  at 
etc . 

Note. that  the  four-signal  test  is  computed  from  10 
beams;  4  beams  of  the  input,  U»\  ,  and  6  of  the  array  response, 
In  general,  the  KA-s  ignal  test  requires  M  beams  of  input 
data  (U)  ,  and  m(m-i)/2  beams  on  the  array  response,  for  a 


T 

Zl 


• .  v;.;  ■ 


total  of  iVl^I/i  r| yZ,  beams  to  compute  the  least-squares  error 


at  any  point  in  the  2M~dimensional  space.  But  the  beams  on 


the  array  response  are  computed  from  the  same  complex  trig¬ 


onometric  terms  that  are  required  for  the  M  beams  of  the  input 
data.  So  the  M-signal  test  requires  evaluation  of  zM  N  sine 
and  cosine  terms  to  compute  the  error  at  any  point  (  N  is  the 


number  of  sensors  in  the  array) .  Thus  the  number  of  trigo¬ 


nometric  terms  requiring  computation  increases  linearly  with 

M. 


It  must  be  noted  that  a  multiple  signal  test  is  not 


everywhere  well-behaved,  but  has  a  singularity.  For  example, 


in  the  case  of  the  2-signal  test,  if 


so  that 


4— k 

Q— P 


E  — i  _  _ 

dzz.  is  undefined.  The  value  it  will  take  on  at  1-k 
depends  on  the  direction  from  which  14.  Though  th 


Though  this 


can,  of  course,  be  shown  analytically,  it  is  a  bit  tedious  for 


repetition  here.  The  contoured  map  of  an  example  (figure  ) 


displays  this  characteristic  graphically.  The  contoured 


function  is  the  2-signal  test 


A  P  +  BCT 


-  ,  , 

^  *  r o  ' 

(a  f-  j/.j^  Q. 

total  of  m(m  *i)/£  beams  to  compute  the  least-squares  error 
at  any  point  in  the  &M-dimensional  space.  But  the  beams  on 
the  array  response  are  computed  from  the  same  complex  trig¬ 
onometric  terms  that  are  required  for  the  M  beams  of  the  input 
data.  So  the  M-s ignal  test  requires  evaluation  of  2MN  s  ine 
and  cosine  terms  to  compute  the  error  at  any  point  (  N  is  the 
number  of  sensors  in  the  array) .  Thus  the  number  of  trigo¬ 
nometric  terms  requiring  computation  increases  linearly  with 

M. 

It  must  be  noted  that  a  multiple  signal  test  is  not 
everywhere  well-behaved,  but  has  a  singularity.  For  example, 
in  the  case  of  the  2-signal  test,  if 


so  that 


k 

a— p 


E  — 1  _  _ 

is  undefined.  The  value  it  will  take  on  at 
depends  on  the  direction  from  which  u.  Though  th 


Though  this 


can,  of  course,  be  shown  analytically,  it  is  a  bit  tedious  for 
repetition  here.  The  contoured  map  of  an  example  (figure  ) 
displays  this  characteristic  graphically.  The  contoured 
function  is  the  2-signal  test 


A  P*+  BCT 


with  !<,  held  fixed  as  varies  over  the  plane.  Note 

that  the  contour  lines  all  run  together  at  =  K, 

'‘i*,  may  range  arbitrarily  close  to  1<  but  must  not 
take  on  that  value  exactly.  The  data  in  this  figure  consist 
of  2  closely  spaced  signals.  The  fixed  vector,  ,  was 

set  at  the  peak  of  their  merged  main  lobes. 

One  might  dismiss  this  singularity  from  practical  con¬ 
sideration  since  signals  of  identical  speed  and  bearing  are 
indistinguishable  by  array  methods.  The  test  for  2  signals  at 
the  same  wavenumber  location  is  thus  unnecessary  anyway.  But 
if  the  2-signal  test,  say,  is  applied  to  data  composed  of  only 
1  signal,  must  not  both  the  probe  vectors  approach  the  same 
point,  i.e.,  the  wavenumber  location  of  the  input  signal,  in 
order  to  merge  and  reduce  the  function  to  the  l-signal  test? 

We  have  seen  that  when  the  data,  ,  consist  of  the  same 

number  of  signals  as  that  for  which  one  is  testing,  the 
test  performs  as  expected:  the  error  is.  minimized  at  the 
wavenumber  location  of  those  input  signals,  and  the  signals 
are  recovered  undistorted.  Suppose,  though,  that  the  2-signal 
test,  say,  is  applied  to  data  consisting  of  just  plane  wave. 


lt  =  Fe1^ 


T  -  It -A 


i  k-i^ 


-  E> 


»*.- 


Let 


iii  Ntitejac  I j  ,  I  iiii  j  (Jfi  tit 


m  the  error  expression 


We  have  to  maximize 


=  ttz:ixi 


A  P*  +  BQ*  . 

ps  p  ,  / 

'  r  N  Ce  ,arvj  ( 

k  qoes  to 


ep  ^s| 

goes  to  ,  then 

P  *  F,  Q  -  FE* 


A  =(P-QE)/(|-E*E) 

becomes 

A  ■*.(  F  -  F  E*E)/(l  -  E*E)  =  F 

and 

6  -(a-PE^I-E^E)  '(FE^-FEVtl-^E) 

-  o 

and 

*■  -  irt  uru,  -  (AP"+  BQ*) 

=  ^^e^KFe^)  -(F"F  +  o) 

=  jr|:F'F  F*F  ■  o 

When  k  goes  to  Ifc  tlvs  error  is  minimized,  the 
signal,  F  ,  is  recovered  undistorted,  and  the  hypothesized 
second  signal  vanishes.  This  solution  is  invariant  though 


^  be  permitted  to  range  over  the  entire  k -plane, 
excepting  the  point  U  .  Thus  the  2-signal  test  does 
not  reduce  to  the  ordinary  f ~ k.  spectrum  in  the  presence 
of  a  single  plane  wave,  and  'tc  is  not  required  to  go  to 
l<,  nor  would  the  gradient  of  £  with  respect  to 
lead  to  l<  (if  one  were  using  a  steepest  descent  technique 
to  minimize  £  ) . 


Numerical  Solution  of  the  Multiple  Signal.  Test 


One  might  propose  to  carry  out  the  numerical  solution 
of  a  multiple  signal  test  by  a  straightforward  search  of 
the  entire  wavenumber  space  of  interest,  as  is  done  in  the 
computation  of  the  conventional  4*  “  spectrum.  But  the 
multiple  signal  test  may  be  used  in  more  practical  fashion, 
with  greater  efficiency,  as  a  follow-up  to  the  ordinary 
4  ”  U  spectrum.  Since  a  high-resolution  array  process 
by  design  is  intended  to  separate  signals  otherwise  unre- 
solvable,  there  is  sound  justification  to  limit  its  use  to 
the  vicinity  of  signals  tentatively  identified  beforehand 
by  less  powerful  but  faster  techniques.  This  is  an  advan¬ 
tageous  circumstance,  since  an  M-signal  test  is  a  function 
of  3M  dimensions  of  wavenumber  and  would  otherwise  prove 
computationally  less  efficient.  Applying  the  X  -signal 
test  to  the  highest  peak  of  an  ordinary  4  ~  i<  spectrum, 
then,  one  hypothesizes  the  presence  of  2  plane  waves  which 
appear  as  1  only  because  of  their  proximity.  By  the  hypothe¬ 
sis  the  spectral  peak  lies  within  the  area  of  the  main  lobe 
of  either  signal  and  thus  ^  may  be  minimized  directly  by 


the  method  of  steepest  descent.  This  is  the  procedure  used 
here . 

Since,  as  has  been  shown  earlier, 

4  s  k 

is  prohibited,  the  descent  cannot  begin  from  any  one  single 
point  in  the  k  -plane,  as,  for  example,  the  peak  under  con¬ 
sideration.  But  any  pair  of  points  in  that  vicinity  is 
suitable;  all  lead  to  the  same  solution.  A  convenient  pair 
are  (1)  the  peak,  and  (2)  the  adjacent  minimum  of  with 
respect  to,  say,  when  is  fixed  at  the  peak  as  in  the 

previously  discussed  figure  .  The  gradient  of  ^  is  com¬ 
puted  at  this  pair  and  itself  then  recomputed  at  a  new 

location  down  the  gradient.  The  length  of  this  first  step 
in  the  descent  is  some  fraction  of  the  width  of  the  array- 
response  main-lobe,  thus  chosen  to  ensure  that  the  process 
does  not  jump  from  the  vicinity  of  the  solution  into  the 
range  of  an  adjacent  relative  minimum.  The  gradient  is  new¬ 
ly  computed  at  this  second  location;  another  somewhat  smaller 
step  is  taken  down  the  gradient;  the  gradient  is  once  more 
computed,  now  at  this  third  location,  and  so  forth  in  suc¬ 
cessively  smaller  steps  until  the  point  is  reached  in  that 
4-d.imensional  space  at  which  the  gradient  goes  to  zero. 


& 


Some  examples  with  synthetic  data  of  the  LASA  LP  array 


follow . 


w 


The  north-south  half-width  of  the  array-response  main- 
lobe  is  about  0.0056  cycle/km.  (figure  )  In  the  first 
example  the  input,  Ua  ,  consists  of  a  signal  at  0.0002 
cycle/km  north,  and  one  at  0.0002  south.  Thus  the  half- 
width  of  the  main  lobe  is  more  than  an  order  of  magnitude 
greater  than  the  distance  between  the  signals  in  the  k.- 
plane.  The  signals  are  aqual,  of  unit  size,  their  phase  is 
equal  (at  the  center  of  the  array),  and  no  noise  is  present. 
The  ordinary  ■f**  k  spectrum,  showing  the  merged  signals 
with  resultant  solitary  main  lobe,  and  looking  precisely 
like  the  array  response,  is  given  in  figure  .  With  one 

vector  fixed  at  the  peak  of  this  main-lobe  while  the  other 
ranges  the  k  -plane,  €  appears  in  the  contoured  plot 

of  figure  .  From  this  pair  of  points,  i.e.,  the  peak  in 

figure  and  the  one  in  figure  ,  the  descent  is  begun. 

Its  progress  and  the  final  result  are  shown  in  the  computer 
bulletin  of  figure  .  This  best  fit  precisely  recovers 
the  2  signals:  size,  phase,  and  the  wavenumber  location  to 
within  less  than  0.000005  cycles/km.  The  wavenumber  distance 
between  these  2  signals  is  only  a  degree  or  two  of  azimuth 
for  Rayleigh  waves. 

In  the  second  example,  presented  in  similar  format 


(figures  and  ),  the  synthetic  input  consists  of  a 
pair  of  signals  in  incoherent  noise. 

One  signal,  of  unit  size  and  zero  phase,  is  located  at 
0.002  cycle/km  south.  Random  numbers  added  to  it  reduce 
the  signal-to-noise  ratio  to  1.  Finally,  a  very  large  signal, 
100  times  the  size  of  the  first,  is  located  at  0.002  cycle/km 
north.  It  is  opposite  in  phase  to  the  first  one  and  thus  it 
interferes  destructively  with  the  small  signal. 

The  location  of  the  2  signals,  superposed  on  the  array 
response,  is  shown  in  figure  .  The  distance  between  them 
is  0.7  the  main-lobe  half-width.  The  arrow  indicates  the 
displacement  of  the  smaller  signal  as  recovered  by  the  2-signal 
test.  The  small  signal  alone .  in  the  presence  of  this  same 
noise  sample,  emerges  with  the  same  displacement  (in  the  ordin¬ 
ary  *£**  k  spectrum)  .  The  steps  in  the  descent  to  the  solu¬ 
tion,  are  presented  in  figure  .  The  distortion  of  amplitude 
and  phase  of  the  large  signal  as  recovered  is  about  1  percent; 
that  of  the  small  signal,  less  than  5  percent.  The  incoher¬ 
ent  noise,  of  course,  is  the  source  of  such  distortion  as  is 
present.  The  modeling  process,  being  linear,  separates  plane 
waves  with  fidelity,  as  demonstrated  in  the  first  example. 
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